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ABSTRACT 

Using CARMA, we imaged the 87 GHz SiO v=0 J=2-l line toward Orion-KL with 0.45" angular 
resolution. The maps indicate that radio source I drives a bipolar outflow into the surrou nding 
molecular cloud along a NE-SW axis, in agreement with the model of Green hill et al.l (|2QQ4l ). The 
extended high velocity outflow from Orion-KL appears to be a continuation of this compact outflow. 
High velocity gas extends farthest along a NW-SE axis, suggesting that the outflow direction changes 
on time scales of a few hundred years. 

Subject headings: ISM: individual(Orion-KL) — ISM: jets and outflows — masers — stars: formation 



1. INTRODUCTION 

Bipolar outflows are ubiquitous from young low-mass 
stars, but are difficult to observe toward high-mass stars, 
which typically form in crowded, physically complex re- 
gions that are quickly disrupted by ionization fronts and 
stellar winds. 

The nearest region of massive star formation is 
the Kleinmann-L ow Nebula in Orion, at a distance 
of ab out 400 pc ([Sandstrom et all 120071 : iMenten etlll 
| 2007f). Measuremen ts of H2O maser proper motions led 
iGenzel et al.l (|1981[ ) to suggest that two distinct out- 
flows originate from this region - a low velocity (18 
km s _1 ) outflow along a NE-SW axis, and a high ve- 
locity (30-100 km s _1 ) flow extending roughly NW- 
SE. Both outflows were inferred to originate within 
a few arcseconds of radio source I, a young star 
with a luminosity of 10 4 to 10 5 L ([Gezari et al.l 
1998). The high velocity o utflow also manifests itself 
as lobes of shock-excited H2 (Be ckwith et al.| [l978) and 
as broad, w eakly bipolar line wings i n CO and other 
mole c ules jZ uckerman e t al.l fl976b Kwan & Sc ovilld 
119761 : lErickson et al.lll982l : IChernin fc Wri ght 1993). 

Source I is associated with a cluster of SiO masers. 
From fits to 2" resolution BIMA data, Plam beck et al.l 
(1990) found that the J=2-l v=l SiO masers were clus- 
tered al ong two arcs offset ~ 0.08" NW and SE of 
source I. IPlambeck et al.l ([19901 ) reproduced this pattern 
with a model of maser emission from a rotating, expand- 
ing disk, tilted at 45° to the plane of the sky. The axis 
of the model disk was projected at PA 145°, suggesting 
that the high velocity outflow emerged along its poles. 

Later observations showed that the J =2 -1 SiO line in 
the v =0 vibrational level also was masing ([Wright et al.l 
119951 ) in a 1000 AU long hourglass-shaped region cen- 
tered on source I. Because the hourglass was elongated 
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NE-SW, it was natural to associate it with the outer 
regions of the model v=l maser disk. SiO linewidths 
were ~ 30 km s _1 across the ent ire hourglass, incon- 
sistent with Keplerian rotation, so Wrigh tet al.l (p. 995) 
suggested that the emission traced the turbulent bound- 
ary layer between an underlying, unseen, disk and the 
high velocity outflow. 

More recent observations have cast doubt on this 
picture. VLB A observations of the 43 GHz v=l 
SiO masers ([Greenhill et al.lll998l : iDoeleman etTaO '1999) 
showed that the brightest maser spots were clustered in 
4 g roups, rather than th an in a ring as predicted by 
the IPlambeck et aT] (p.990) model. A brid ge of emission 
conne cting the two southern clusters led I Greenhill et al.l 
(2004) to propose that the v=l masers originate in a 
nearly edge-on disk rotating about a NE-SW axis, per- 
pendicular to the previously hypothesized disk, and that 
the v=0 SiO emission traces the base of the 18 km s _1 
outflow. This model did not attempt to explain the high 
velocity outflow. 

Proper moti on m easurements of source I 

([Rodriguez et al.l 120051 : iGomez et al.l 120081 ) also pose 
problems for the original model. These data indicate 
that source I is moving to the SE at 0.007"/year, 
plowing through the molecular cloud at ~ 14 km s _1 . 
In that case a 1000 AU diameter disk would rapidly be 
stripped away by the ram pressure of the ambient gas 
unless source I were implausibly massive. 

In this paper we report 0.45" (180 AU) resolution im- 
ages of the SiO v=0 J=2-l line obtained with CARMA. 
The new maps provide the best evidence to date that 
the v=0 SiO emission or iginates in a NE-SW outflow 
from source I, as in the Greenhi ll et al.l (|2004f ) model. 
The maps also suggest that the high velocity outflow is 
a continuation of this compact flow. 

2. OBSERVATIONS 

Observations were made with CARMA in the A, B, 
and C arrays, providing projected antenna spacings rang- 
ing from 5 to 545 kA. The total integration time on 
Orion-KL was 3.4 hrs in the A-array (2009 Jan), 3.7 hrs 
in the B-array (2008 Feb), and 5 hrs in the C-array (2009 
May). 

The correlator was configured for simultaneous obser- 
vations of the J=2-l SiO transitions in both the v=0 
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(86.847 GHz) and v=l (86.243 GHz) vibrational levels. 
For each transition the velocity coverage was 104 km s _1 
and the resolution was 3.4 km s _1 after Hanning smooth- 
ing. Observations of quasars were used to calibrate phase 
and amplitude ripples across the I.F. passband. After 
these were removed, the visibility data were selfcalibrated 
using a 28 km s _1 wide channel that included all the 
bright v=l maser features. Only baselines shorter than 
250 kA were used in the selfcalibration because the maser 
is slightly resolved on the longest baselines. The integra- 
tion time per data record was 4 sec in the A-array and 
10 sec in the B and C arrays, allowing near-complete re- 
moval of atmospheric phase fluctuations; phase residuals 
on the maser were typically ±1° on all baselines after 
selfcalibration. 

Maps of the v=0 SiO line were made with the Miriad 
data reduction package. Uniformly weighting the visi- 
bility data yielded an 0.54" x 0.40" synthesized beam at 
PA 52°. Because the data were selfcalibrated using the 
v=l SiO masers, the maps are centered on the intensity- 
weighted mean maser position, which does not coincide 
precisely with radio source I. An 86 GHz continuum 
map made with A-array data from 2009 Jan showed that 
source I was centered 0.01 ±0.01" W and 0.03±0.01" S of 
the mean v=l maser position. We adjusted the central 
coordinate of the maps to place source I at 05 /l 35 m 14!515, 
— 5°22 / 30"57, as d erived from VLA proper motion data 
(jGomez et al.ll2008[ ). 

3. RESULTS 

Figure [T] presents a series of channel maps showing the 
v=0 SiO emission in an 8" x 8" box centered on source I. 
The lowest contour on these images is 75 K and the peak 
is 3800 K. 

Although the maps in Figure [T] a re con sistent 
with those published by I Wright et al.l ((1995) and 
iChandler de Pred (|1995f ). the higher quality CARMA 
images clearly indicate that v=0 SiO emission origi- 
nates in a NE-SW bipola r outflow from source I, as in 
the iGreenhill et al.l ([2004) model. The central channels 
show the limb-brightened edges of two cones centered 
on source I; emission in the line wings is brightest in- 
terior to these cones. Although the outflow appears to 
lie nearly in the plane of the sky, there is a measurable 
red-blue asymmetry: redshifted gas at Vlsr= 22 km s _1 
appears as a narrow jet to the SW, while blueshifted gas 
at Vlsr= — 9 km s _1 is offset to the NE. This asymmetry 
is inconsistent with the old model in which the emission 
originates in the turbulent boundary layer of a disk. 

It is clear that the SiO outflow corres ponds to the 
18 km s _1 H2O maser outflow identified bv lGenzel et al.l 
(1981): toward source I the SiO line has a fullwidth at 
zero intensity of ~ 36 km s _1 , and both SiO and H2O 
masers extend along a NE-SW axis. It is likely, however, 
that this material was accelerated by a much faster wind. 
Source I is traveling to the SE at 14 km s _1 in the plane 
of the sky (jGomez et al.l l2QQ8h . so gas moving outward 
at only 18 km s _1 would appear to be swept back into 
an arc with arms trailing to the NW. 

The high signal to noise level in the maps, exceeding 
100:1 on the bright maser spots, makes it feasible to en- 
hance the resolution of the innermost region of the out- 
flow. One can fit the position of an isol ated maser featur e 
to an accuracy of ~ 0.5 #fwhm/SNR (Rei d et aLl ll988). 
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Figure 1. Channel maps, 5.06 km s _1 wide, of the SiO v=0 J=2- 
1 line toward Orion-KL. Each box is 8 /x on a side and is labeled 
with the center LSR velocity. A cross marks the location of radio 
source I. The contour levels are ± 75, 150, 250, 400, 600 K, then 
900 to 3400 K in 500 K steps; 745 K = 1 Jy/beam. The noise level 
is 12 K, increasing to 19 K in the center channels due to limited 
dynamic range. 

where #fwhm is the apparent FWHM of the source and 
SNR is the signal-to-noise ratio. Here the maser spots 
overlap heavily, however, so fitting their positions indi- 
vidually is difficult. Instead, we CLEANed the maps in 
the usual way, then convolved the CLEAN components 
with a 0.25" FWHM Gaussian restoring beam, about 
half the width of the synthesized beam. Figure [2] shows 
the resulting SiO v=0 image at Vlsr= 5 km s _1 in the 
2" x 2" central region. Red contours show the 229 GHz 
continuum emission from source I, which is centered in 
the waist of the SiO hourglass; SiO masers in the v=l 
vibrational level are clustered along the edges of the con- 
tinuum source. 



Outflow from Orion Source I 
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Figure 3. Velocity moment map of the v=0 J =2-1 SiO line de- 
rived from \" resolution channel maps covering the velocity range 
—30 <Vlsr< 47 km s _1 . Pixels were averaged into the mo- 
ment map only if they exceeded 3x the rms noise, which ranged 
from 4-8 K across the channels. Black contours show the Vlsr= 
6.7 km s — 1 channel map from Figure 1, with the same contour in- 
tervals. A white box marks the location of the Becklin-Neugebauer 
Object; radio continuum from BN corrupts the velocity moment 
map in this direction. White arrows show proper moti ons of the 
HH o bjects 152-228, 152-229, and 135-233 over 150 years (|Doi et al.l 

Eqq2). 



Figure 2. (top) SiO v=0 map in a 10 km s 1 wide channel 
centered at Vlsr= 5 km s _1 , deconvolved to 0.25" resolution as 
described in the text. The contour interval is 780 K (0.3 Jy/beam), 
the rms noise is 80 K, and the peak intensity is 6250 K. Red con- 
tours show the 229 GHz radio continuum from source I, measured 
at CARMA with an 0.15 /x synthesized beam (Plambeck et al., in 
preparation), (bottom) Position- velocity cuts along lines (a) and 
(b) generated from 3.4 km s _1 resolution channel maps, also de- 
convolved to 0.25" resolution. The contour spacing is 780 K; the 
rms noise of the channel maps ranges from 50-140 K. Green el- 
lipses show velocities expected for a 0.5" diameter ring expanding 
radially at 10.5 km s _1 in panel (a), 9 km s _1 in panel (b). 



The brightest v=0 SiO masers lie along two bars offset 
~ 0.5" NE and SW of the continuum source. The ends of 
the bars are limb-brightened, suggesting that the masers 
originate in two annuli. Position-velocity cuts through 
the bars, shown in panels (a) and (b), hint that these 
annuli are expanding radially at ~ 10 km s _1 , as modeled 
by the green ellipses. 

What about the high velocity outflow? In order to ob- 
tain higher sensitivity for extended emission, we tapered 
the weighting of the visibility data to produce a set of 
channel maps with a 1.13" x 0.93" synthesized beam. 
The rms noise is 4 K, increasing to as much as 8 K in 
channels with strong maser emission due to limited dy- 
namic range. Figure [3] shows the full extent of the SiO 
emission detected in these maps, color coded to indicate 
the intensity-weighted LSR velocity. Brightness temper- 
atures outside the central 8" x 5" hourglass are < 100 K, 
hence this emission is likely to be the rmal. As in previ- 
ous m aps of the high velocity outflow (jChernin fc Wrightl 
1996), the strongest redshifted emission is E of source I, 
while blueshifted gas is offset to the NW. 

Figure [3] suggests that the extended high velocity out- 
flow is simply a continuation of the compact outflow. 
For example, ~ 3" NE of source I the compact outflow 
twists to the E and points directly toward some of the 



most redshifted gas in the extended flow. The opening 
angle appears to broaden abruptly at this point, per- 
haps as the outflow breaks out of dense gas surrounding 
source I. Two HH objects, 152-228 and 152-229, lie at 
the tip of the flow; proper motion measurements indicate 
that they are mo ving away from source I at 35-50 km s _1 
(|Doi et al.ll2002h . 
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Figure 4. SiO v=0 channel maps showing a velocity gradient 
in the redshifted lobe of the high velocity outflow. The contour 
interval is 12 K; the rms noise is 4 K. The Vlsr of each channel 
is indicated. Gray shading indicates pixels for which there is 3cr 
emission in at least one channel in the velocity range 29 <Vlsr< 
52 km s — 1 . Blue stars mark the positions of source I and BN. 
Red arrows show proper motions of the HH objects 152-228 and 
152-229. 



Figure [3] suggests that SiO velocities increase away 
from source I. This is particularly evident in the red- 
shifted £ as 5 as shown in Figure [4j note the absence of 
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high velocity gas within 4" of the star. There are several 
possible explanations for these velocity gradients. 

First, SiO emission may originate in dense clumps that 
are entrained in a fast stellar wind. The wind continues 
to accelerate the clumps as they move outward. Evidence 
for such acce leration is relatively common in outflows 
(|StaHerill994h . 

Second, it's possible that the SiO emitting gas was 
ejected by an explosive event that flung out material with 
a range of speeds. The fastest-moving ejecta now are 
farthest from the center. The system of shock-excited H2 
fingers in Orion appears to have been created by such a n 
explosive event ([Allen fc Burtonl 119931 : iDoi et al.ll2QQ2h . 
perhaps as a swept-up shell was fragmented by a faster- 
moving wind ([Stone et al.lll995l ). 

Third, it's possible that the extended gas provides 
a fossil record of a precessing outflow. Currently the 
outflow axis is NE-SW, with the (redshifted) SW lobe 
tipped slightly into the plane of the sky. Approximately 
10" from the star, the outflow axis appears to be E-W, 
with the E lobe tipped into the plane of the sky. If the 
true outflow velocity is 100 km s _1 , this > 45° change in 
direction took place in 200 years. Precession occurs be- 
cause of tidal interactions in noncoplanar binary systems; 
periods of a few thousand years are typ ical of systems 
with binary separations of tens of AU (|Ter quern et al.l 
1999). A secondary in an eccentric orbit can change both 
the outflow direction and the mass lo ss rate periodically, 
as in t he model of CepA presented by I Cunningham et al I 
(2009L 



Zapata et al.l (|2009l ) report that much of the emission 



from the high velocity CO line wings in Orion originates 
from filamentary structures similar to th e H2 fingers, 
with linear velocity gradients along them. Za pata et al.l 
(2009) argue that the high velocity flow is not a clas- 
sical stellar outflow, but instead was produced by the 
disintegration, about 500 years ago, of a multiple stel- 
lar system that contained source I and BN. While this 
scenario is consistent with the velocity gradients evident 
in Figure [4] and with the vaguely filamentary character 
of the SiO emission away from source I, it implies that 
the apparently continuous transition from the 18 km s _1 
flow to the high velocity outflow seen in Figure [3] is il- 
lusory. We are reluctant to accept that conclusion, and 
suggest that perhaps the disintegration of the multiple 
system triggered an abrupt increase in the velocity or 
mass loss rate of the outflow fro m source I, creatin g the 
finger system as in the model of iStone et al.l ([1995). 

4. CONCLUSIONS 

New maps of the SiO line in the v=0 J=2-l transition 
toward Orion-KL provide some of the clearest evidence 
to date that radio source I, thought to be a massive young 
star, drives a bipolar outflow into the surrounding molec- 
ular cloud along a NE-SW axis. This outflow contains 
multiple water masers and conventionally is referred to 
as the 18 km s _1 flow. The outflow velocity must be 
substantially greater than 18 km s _1 , however, or the 
outflow lobes would trail back to the NW owing to the 
proper motion of source I. Close to source I the SiO v=0 



emission is masing. The strongest masers are clustered 
in two annuli offset ~ 200 AU from source I along the 
central axis of the outflow. 

More extended, higher velocity gas extends along a 
NW-SE axis, almost perpendicular to the 18 km s _1 flow. 
The SiO maps suggest that this weakly bipolar outflow 
is an extension of the 18 km s _1 flow, which may change 
direction on time scales of a few hundred years. 
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Science Foundation under a cooperative agreement, and 
by the CARMA partner universities. 
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